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ABSTRACT 

Intersection with the debris of a large (50-100 km) short-period comet during the 
Upper Palaeohthic provides a satisfactory explanation for the catastrophe of celestial 
origin which has been postulated to have occurred around 12900 BP, and which pre- 
saged a return to ice age conditions of duration ~1300 years. The Taurid Complex 
appears to be the debris of this erstwhile comet; it includes at least 19 of the brightest 
near-Earth objects. Sub-kilometre bodies in meteor streams may present the greatest 
regional impact hazard on timescales of human concern. 
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1 INTRODUCTION 

The sudden onset of the Younger Dryas cooling 12,900 years 
ago was marked by intense wildfires over North America, 
major disruption of human culture, and the rapid extinc- 
tion of 35 genera of North American mammals (Faith & 
Surovell 2009). A thin carbon-rich black layer of this age has 
been identified at many sites across North America (Haynes 

2008) , coincident in age with the Younger Dryas boundary. 
Recently, several geochemical markers at this layer have 

been presented which seem to indicate that a major extrater- 
restrial event was involved in these events (Firestone et al. 
2007; Firestone 2009). The markers include: 

Nanodiamonds, found in abundance at fifteen Younger 
Dryas Boundary (YDB) sites across North America and 
north-west Europe, covering a quarter of the Earth's surface. 
They are mostly found embedded within particles of melted 
plant resins, confined to the thin black layer, none having 
yet been found in other strata dating from 55,000 EC to 
the present. All three diamond allotropes are present in the 
boundary sediments including lonsdaleite (hexagonal nan- 
odiamonds) , which is shock-synthesised and found on Earth 
only in association with ET impacts or inside meteorites 
(Kennett et al. 2009). The nanodiamonds and lonsdaleite 
found in meteorites were already present before atmospheric 
entry (Clarke et al. 1981). 

Soot, micro spherules and magnetic grams found at high 
concentrations, the microspherules having trace element 
abundances comparable to lunar KREEP but not to any 
other terrestrial or observed meteoritic source (Firestone 

2009) . The only previously known co-occurrence of soot, 
nanodiamonds and rapid extinction is at the Cretaceous- 
Tertiary boundary layer. 
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To cause destruction on a continental scale by a single 
impact, the energy required is ~10'' Mt, corresponding to a 
4 km wide comet (Toon et al. 1997). Firestone et al. (2007) 
proposed that such a comet broke up on atmospheric en- 
try and struck the 2-km thick Laurentide ice sheet. They 
speculated that four deep holes in the Great Lakes might 
be craters produced by this event. This short-term catastro- 
phe was followed by a sudden return to ice age conditions 
on Earth which lasted for over a millennium, the cause of 
which is a matter for speculation (loc. cit.; Franzen & Cropp 
2007) . An alternative proposal to explain the data is that the 
Earth encountered a rare swarm of carbonaceous chondrites 
or comets, yielding multiple airbursts and possible surface 
impacts (Kennett, Kennett, West et al. 2009). 

The nature of the event at the Younger Dryas Bound- 
ary remains controversial. For example it has been argued 
that the mammoths were in decline for about a millennium 
before the final extinctions (Gill, Williams, Jackson et al. 
2009) . Nevertheless it seems difficult to account for the geo- 
physical markers - nanodiamonds, exotic sediment composi- 
tion, evidence of a sudden, continent-wide fiery catastrophe 
etc - without the occurrence of some extraordinary extrater- 
restrial event. 

Problems have however been raised from the astronom- 
ical perspective. It is difficult to reconcile a major astro- 
nomical catastrophe of continental dimensions as recently 
as 12.9 ka ago with the known population of near-Earth ob- 
jects (NEOs) currently being revealed by comet and aster- 
oid search programmes such as the LINEAR and Catalina 
surveys. Thus it has been estimated that a 10 megaton im- 
pact may occur anywhere on Earth with a mean recurrence 
time 2000-3000 years (Stuart & Binzel 2004). An impact of 
energy 1000 Mt or more is expected with recurrence time 
~60,000 yr (loc. cit. and Morbidelli, Jedicke, Bottke et al. 
2002). While this would be immensely damaging, the ef- 
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feet would be regional rather than continental or global. On 
conventional reckoning, the impact of a 4 km asteroid as 
recently as 13,000 years ago is a one in a thousand event, 
and if an active comet, the probability is a few in a million 
(Harris 2008). It has also been pointed out that, to ignite 
fires in southern Arizona by thermal radiation from a fire- 
ball over the Great Lakes, the projectile would have had to 
be 100 km across, and all non-microbial life on Earth would 
have been wiped out (Melosh 2009). The geophysical evi- 
dence for an astronomical catastrophe at this boundary has 
been extended to Venezuela (Mahaney et al. 2010), which 
exacerbates that difficulty. 

However, a potential deficiency of impact assessments 
based on current Spaceguard data is that they assume sta- 
tistical completeness. To estimate impact risks at say the 
~10~* per annum level by extrapolating from ~10 yr of 
Spaceguard surveys is analogous to forecasting three years' 
weather by extrapolating one day's observation. A cascade of 
comet disintegration, for example, could be a significant haz- 
ard on such timescales (cf the 1994 impacts on Jupiter fol- 
lowing the disintegration of D/1993 F2 Shoemaker- Levy 9). 

It has also been objected that, if the Younger Dryas 
event was due to an encounter with a swarm of debris 13,000 
years ago, the material would still be visible in the inner so- 
lar system at the present time. Meteoroid data, being largely 
tracers of comet disintegration, do indeed have the poten- 
tial to look into the history of the cometary environment for 
tens of thousands of years, and in this paper we discuss such 
evidence. We first, in Section 2, summarise evidence which 
indicates that a large, short-period comet entered the inner 
planetary system a few 10* yr ago and has since then been 
undergoing a hierarchy of disintegrations, yielding the Tau- 
rid Complex. We show that at least 19 of the largest NEOs 
have orbits significantly close to that of Comet Encke. In 
Section 3 we numerically model the disintegration history 
of the progenitor comet. This reveals that encounters with 
dense swarms of material, sufficient to produce a 12.9 ka 
cosmic event, are indeed reasonable expectations of recent 
Earth history. Section 4 gives a broadbrush description of 
the likely consequences of an encounter with such a swarm. 
There is satisfactory agreement with the main physical fea- 
tures at this boundary. 

The occurrence of a major celestial catastrophe over 
continental dimensions as recently as 12900 years ago, which 
is quite unexpected from Spaceguard data alone, may carry 
implications for assessments of the impact hazard currently 
faced by civilization (cf Napier & Asher 2009). 

2 A LARGE PALAEOLITHIC SHORT-PERIOD 
COMET 

2.1 The zodiacal cloud and the Taurid meteors 

The complex nature of the Arietid- Taurid meteor streams 
has long been recognized, with for example Denning (1928) 
identifying 13 separate sub-streams within the system. The 
Complex encompasses comet 2P/Encke with orbital period 
3.3 yr. Differential orbital precession arising from Jovian 
perturbations yields two main night-time branches, and two 
daytime branches. From the orbital dispersion of the mate- 
rial, Steel & Asher (1996) derive an age 20-30 kyr for the 
system (cf Babadzhanov et al. 1990). 



This broad complex, which comprises at least half the 
mass of the zodiacal cloud, appears not to be in equilibrium 
and to be replenished sporadically (Napier 2001). Observa- 
tions of both visual (Stohl 1986) and radar 'sporadic' me- 
teors (Campbell-Brown 2008) indicate that they are not in 
fact sporadic but are rather concentrated in a small number 
of primary sources, two of which - the helion and and antihe- 
lion sources - appear to derive from Comet 2P/Encke or its 
parent. Wiegert, Vaubaillon & Campbell-Brown (2009a,b), 
modelling this sporadic meteor complex, find that this comet 
dominates the helion and antihelion sources by a factor of 
about ten. They consider that the sporadic population is 
dominated by a small number of comets (especially Encke 
and 55P/Tempel-Tuttle) with high-transfer efficiencies to 
Earth, rather than the comet population as a whole. The 
current Comet Encke is, however, inadequate to supply the 
zodiacal cloud, and it is probably just a recently reactivated 
fragment of the progenitor. 

The Taurid Complex has also been detected in situ 
through the Long Duration Exposure Facility, in the mete- 
oroid mass range 10~^'^ to 10^* gm, on four large-area faces 
over a period of 5.8 years (McBride et al. 1995). These work- 
ers found that the Taurid Complex is a major part of the 
current zodiacal cloud, yielding a significant proportion - say 
a third - of the sporadic meteoroid influx. Zook (2001) found 
that asteroidal grains probably do not contribute more than 
25% of the meteoroid cratering population on LDEF, comet 
dust being the rest. These results appear to strengthen the 
earlier conclusions of Whipple (1967) and others that the 
Taurid meteors are the remnants of an erstwhile exceptional 
body in a short-period. Earth-crossing orbit. 

They are also consistent with the empirical flnding that 
the current zodiacal cloud is overmassive by one or two pow- 
ers of 10 in relation to current replenishing sources (Whipple 
1967; Clube & Napier 1984; Hughes 1996; Napier 2001). The 
mass of the zodiacal cloud is estimated to be ~2-30x 10^^ gm 
(Hughes 1996). If say half of this is due to a single progenitor 
comet, and (conservatively) half its input has disappeared 
due to Poynting-Robertson and coUisional effects, then for 
progenitor comet density 0.5 gmcm^'', this yields an original 
progenitor comet 50 - 100 km in diameter, neglecting the 
contribution from larger bodies within the Complex (Napier 
2009). 

2.2 Near-Earth objects associated with Comet 
Encke 

Comet evolution commonly proceeds via fragmentation into 
smaller bodies, and indeed this may be the prime disinte- 
gration route (Jenniskens 2008; Levasseur-Regourd, Hadam- 
cik, Desvoivres et al. 2009). The disintegration of D/1993 
F2 Shoemaker-Levy 9 was followed by the impact of ~20 
fragments on to Jupiter, with impact energies ~10^ Mt. 
In late 1995, Comet 73P/Schwassmann-Wachmann 3 split 
into three fragments, having been seen as a single object 
on its previous return in 1990. On subsequent returns fresh 
mini-comets were seen, and by 2006 over 150 fragments had 
been detected (Reach, Vaubaillon, Kelley et al. 2009; Ishig- 
uro, Usui, Sarugaku et al. 2009). At least one fragment was 
~300 m across and a dozen or so were in the size range of the 
Tunguska impactor. Comet C/1996 B2 Hyakutake split into 
seven sub-nuclei whose sizes seem to have been in the range 
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Figure 1. Distribution of tlie longitude of perihelion m of the 
19 brightest NEOs {H < 16.5) with (a, e, i) in the range 1.85< 
a <2.7, 0.65< e <1.0 and < i < 14° encompassing the Taurids 
and Comet Encke, but with no restriction on ro. These are shown 
as asterisks. Their mean ra=174° ± 12°, which is 13° from Comet 
Encke's longitude of perihelion m ~161°, marked with a square. 
The inner arc gives the spread of the recognized Taurid meteor 
showers (Stohl 1986). The bottom five asterisks appear to form 
a separate group, which includes the large near-Earth asteroid 
Hephaistos. 
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departures of mean longitude of perifielion from 161 deg 

Figure 2. Outcome of 50,000 trials in which the longitude of 
perihelion ro of 19 objects is randomly distributed and the mean 
departure from ro = 161° — that of Comet Encke - is recorded. 
For comparison the mean absolute difference of 13° is shown, 
derived from the 19 NEOs of Fig. [T] Fragmentation of a large 
comet is inferred to be the source of these orbital correlations, on 
a timcscale '~20-30 kyr. 



of the Tunguska impactor (Desvoivres, Klinger, Levasseur- 
Regourd et al., 2000); Comet C/1999 S4 LINEAR disin- 
tegrated into thousands of fragments. Di Sisto, Fernandez 
& Brunini (2009) consider that on average a single Jupiter 
Family comet (P <20 yr) will undergo a major splitting ev- 
ery < 77 revolutions. There appears to be more than one 
mode of disintegration, breakup at perihelion presumably 
due to tidal stresses being one, and 'unknown' being the 
others. In the case of the Kreutz sungrazers, several hundred 



fragments are known of a total probably running to thou- 
sands, and it appears that disintegration proceeds through 
a cascade of fragmentations (Sekanina & Chodas 2004), the 
fragments themselves continuing to break up. 

Cometary breakup may also involve the creation of dor- 
mant bodies (Kresak & Kresakova, 1987; Jenniskens loc. cit., 
Porubcan, Kornos & Williams 2004; Babadzhanov, Williams 
& Kokhirova 2008a). This seems to have happened, for ex- 
ample, in the case of the Quadrantids (Jenniskens 2004; 
Babadzhanov et al. 2004), and Comet Biela (Jennisken & 
Vaubaillon 2007), while fireball outbursts from the k Cyg- 
nid shower support the hypothesis that the formation of the 
stream and dormant bodies within it proceeded from a cas- 
cade of disruptions of the progenitor comet (Trigo- Rodriguez 
et al. 2009). 

The presence of substreams in the Taurid Complex indi- 
cates that the progenitor comet also took the route of mul- 
tiple disintegration, and searches have revealed a number 
of NEOs associated with the sub-streams of the complex. 
(Clube & Napier 1984, 1990; Asher, Clube & Steel 1993; 
Babadzhanov, Williams & Kokhirova 2008b, Porubcan, Ko- 
rnos & Wilhams 2006). Clube & Napier (1984) pointed to 
orbital resemblances between Comet Encke, Oljato, 1982 
TA and 1984 KB - later numbered and named 6063 Jason 
- and proposed that these bodies, along with the Taurid 
meteors, were the remnants of an erstwhile giant comet. 
These are amongst the largest NEOs. The database has ex- 
panded enormously since then, and can be further examined 
to see whether the hypothesis continues to hold up. The 
data chosen for this purpose were the NEAs listed in the 
Minor Planet Center to 2nd April 2008, comprising 2250 
Amor (Mars-crossing) asteroids and 2622 ApoUos (Earth- 
crossing). 

Fig. [T] shows the result of extracting the 19 brightest 
near-Earth asteroids {H < 16.5) in the range 1.83< a < 
2.7 AU, 0.63< e < 1.0 and 0° < i < 14° broadly encompass- 
ing the elements of Comet Encke, but with no restriction 
on the longitude of perihelion w, which is --^161° for the 
comet. In the absence of a generic relation between these 
NEAs, their tu-distribution is expected to be random, with 
the longitudes of perihelion precessing at somewhat differ- 
ent rates. However it can be seen that this is not so, the 
average value of w being ~179±14°, i.e. 13° from Comet 
Encke, as against an expectation of ~90° for a random dis- 
tribution (Fig. [H). AU but one of these NEOs lie in the 
120° range 133< uj < 251° encompassing Comet Encke. 
The significance of this result was tested by running trials 
in which the longitude of perihelion of 19 objects was ran- 
domly distributed, and their mean departure from 161° was 
tested. The outcome of 50,000 such simulations is illustrated 
in Fig. (2] Clearly, the configuration shown in Fig.[T]has neg- 
ligible probability of being due to chance. 

Five of the 19 NEAs in Table [T] are tightly bunched 
around the current w — 237° of 2212 Hephaistos, their mean 
longitude being 238±7.6°. They appear to form a separate 
group (Clube & Napier 1984; Steel & Asher 1994). If these 
and Hephaistos are excluded, the remaining NEAs have 
<-cu >= 140. 8± 10.1° as against say ~150.5° for the 15 Tau- 
rid substreams identified by Porubcan, Kornos & Williams 
(2006), or 136.2° for the 10 Taurid Complex asteroids which 
they identified. The longitudes of perihelion of the Taurid 
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Table 1. Bright, high eccentricity NEAs selected as described 
in the text. The longitudes of perihelion tend to cluster in two 
groups corresponding to the Taurid meteor complex and 2212 
Hephaistos respectively. 



name 


u 

rl 


1a7 


i 


e 


CL 


1996 FR3 


16.40 


63.9 


8.1 


0.795 


2.167 


Poseidon 


15.50 


123.8 


11.9 


0.679 


1.836 


1982 TA 


14.60 


129.4 


12.2 


0.772 


2.302 


Lugh 


16.40 


133.0 


3.9 


0.698 


2.567 


Zeus 




LOO.O 


114 
± 1 


U.UO'rt 




5025 P-L 


15.70 


139.0 


3.8 


0.746 


2.523 


1998 SS49 


15.80 


143.9 


10.8 


0.640 


1.923 


Jason 


15.30 


146.5 


4.9 


0.763 


2.220 


1999 UM3 


16.30 


150.3 


10.7 


0.669 


2.383 


Cuno 


14.40 


171.0 


6.8 


0.636 


1.981 


Oljato 


15.20 


172.8 


2.5 


0.713 


2.172 


Heracles 


14.00 


177.2 


9.1 


0.773 


1.833 


1994 AH2 


16.30 


189.2 


9.6 


0.709 


2.533 


2000 WK63 


16.20 


205.0 


10.4 


0.758 


2.438 


Hephaistos 


13.90 


236.8 


11.7 


0.834 


2.167 


1990 TGI 


14.70 


238.8 


8.7 


0.680 


2.440 


1990 SM 


16.10 


243.9 


11.6 


0.766 


2.101 


1999 RD32 


16.40 


250.0 


6.8 


0.771 


2.638 


Mithra 


15.60 


251.1 


3.0 


0.661 


2.204 



meteors generally lie in the range 100° < < 190° (e.g. 
Stohl 1986). 

A possible objection to the idea of a common origin for 
these bodies is that, where the surface reflectance properties 
have been measured, they are more akin to common main 
belt asteroids than the D-type, carbonaceous surfaces associ- 
ated with comet nuclei and the more volatile rich, 'primitive' 
regions of the main asteroid belt. While there is spectral di- 
versity amongst the nuclei of comets (Luu 1993), not enough 
is known about the spectral characteristics of comet nuclei 
to permit a confident diagnosis about the provenance of any 
particular small body. Nearly all the NEAs in Table [T] have 
perihelia <1 AU, and thermal histories quite different from 
those of most comets. Sublimation of volatiles at high tem- 
peratures from the surface layers may leave behind mantles 
of rubble relatively depleted in dark organics. A possible 
example of this is the 5 km object Phaethon, which has 
the reflectance spectrum of a main belt asteroid, but whose 
associated Geminid meteors give them an unambiguously 
cometary origin (Kasuga 2009; Borovigka 2007). Ten of the 
Taurid Complex asteroids listed in Table [l] also appear in a 
group of 17 NEAs obtained by Babadzhanov (2001) using 
the reduced D-criterion of Asher et al. (1993), when orbital 
similarity was tested against that of 2P/Encke. The list of 
Babadzhanov (2001) was not restricted by magnitude. All 
the NEAs in his list had meteor showers associated with 
them, supporting a cometary origin for these bodies. 

The alignment shown in Fig. [l] is rapidly lost as one 
moves to fainter magnitudes. This might imply that, at the 
sub-kilometre level, there is a dearth of bodies produced 
by the disintegration of the progenitor comet, or that its 
signal is lost in the noise of asteroids from the main belt. 
Detection of such bodies is likely to be difficult, although 
Porubcan et al. (2006) and Babadzhanov et al. (2008b) have 
found dynamical evidence for nine small asteroids, as well 
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Figure 3. A 1,000-particlc meteoroid trail viewed from an eclip- 
tic pole. The trail was generated by fragmontating a comet in 
an Encko-like orbit at perihelion, the fragments being ejected in 
random directions at speeds up to 20 ms~^. The Earth would 
traverse this trail in at most a few hours. 

as 15 substreams, within the Taurid meteor complex and 
associated with observable meteor showers and fireballs. 

In the course of its disintegration, one expects the orig- 
inal comet to have given rise to concentrations of material 
along its orbital track. This disintegration may be contin- 
uing to the present day, but given the evidence that the 
original comet was exceptionally massive, there may have 
been epochs of strongly enhanced fireball and Tunguska-like 
activity during its most active phase. Swarms of particles 
trapped in resonant orbits would extend the lifetime of some 
of these concentrations, as predicted by Asher and Clube 
(1993) and observed by Dubietis & Arh (2007). Porubcan 
et al. (2006) found that two of the filaments originated 4000 
- 5000 years ago, while one originated over 5000 years ago. 
The Complex is likely to be much older than this, perhaps 
0.2-0.3 Myr (Steel & Asher 1996). This great age, deduced 
from the dispersion of the material, is consistent with its ex- 
ceptional mass. Since the putative 12,900 BP event of North 
America fits comfortably within the timescale of evolution 
of the progenitor comet, we may reasonably ask whether one 
can construct a comet disintegration model to account for 
the event. 



3 EVOLUTION OF A LARGE 
DISINTEGRATING COMET 

Generally, the separation speeds of cometary fragments are 
up to ~5 ms~^ (Kresak 1993). However the escape veloc- 
ity from the surface of a 100 km comet is ~15 ms~^. As 
the bodies cascade down through smaller sizes, disruption 
speeds will presumably decline. Here we consider the case 
where a comet fragment in an Encke-like orbit disrupts at 
perihelion. The initial orbit is given an inclination of 6° . The 
fragments are assumed to move away in random directions 
from the comet. 

Figure |3] shows the projection, looking from an eclip- 
tic pole, of 1,000 particles ejected at perihelion with ran- 
dom speeds up to Vm ~ 20 ms~^, during their first passage 
through a heliocentric sphere of radius 1 AU. Depending on 
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the encounter geometry, the Earth traverses this trail in at 
most a few hours. For comparison Kresak (1993) found that 
the duration of rare, intense meteor storms is about an hour. 
Tlie cross-sectional area of this trail is ~2000 times that of 
the Earth, but its length increases linearly with time. 

The Opik encounter formula (Opik 1976) leads to an 
expectation that, for a single object in an orbit like that of 
Comet Encke, the mean interval between collisions with the 
Earth is ~2xl0* yr. For an encounter with a single mete- 
oroid trail with lO'' times the area of the Earth, the interval 
between such encounters becomes ~2 x 10^ yr, comparable 
with the 20-30 kyr fragmentation history of the Taurid Com- 
plex progenitor. Thus if at any given time - during its major 
disintegration history - such a trail is likely to be present, 
one or more such encounters is a reasonable expectation. If 
the meteoroid trail comes from the perihelion disintegration 
of say a 5 X 10^^ gm fragment of the original comet, then 
the Earth would encounter say 1-5x10^'* gm of this mate- 
rial, energetically equivalent to the impact of 2,000-10,000 
Tunguska objects over about an hour. 

The fragments from a disrupting comet tend to have 
size distributions oc d~'' with index q not very different from 
those expected from self-similar cascades {q = 3.5), imply- 
ing that the mass tends to be concentrated in a few largest 
bodies (e.g. Ishiguro et al. 2009). The probability of a dam- 
aging encounter was investigated by adopting a self-similar 
fragmentation model, in which half the mass of a disinte- 
grating fragment was assumed to be lost as dust, and the 
mass of the largest sub-fragment was taken as 0.3 times the 
initial mass. This yielded five sub-fragments, each of which 
disintegrated in self similar manner. Assume that the limit- 
ing mass for a meteoroid trail - generated at perihelion and 
intercepted at Earth - to cause a damaging bombardment is 
rric = lO^'^ gm. Each cometary fragment poses a hazard un- 
til it disintegrates into sub-fragments and so on, hence the 
totality of trails created during a cascading disintegration 
must be computed. With this model we find that a comet 
whose initial mass is ~10'^'' gm will, in the course of its disin- 
tegration, yield ~1,000 trails of mass at least 10^^ gm. Then 
by the above argument, if the active lifetime of the comet 
is say 20,000 years, the swarm must remain as a coherent 
entity for ^^20 yr for an expectation of one such encounter. 

This lifetime requirement is modest, the more so since 
the orbital period of the Taurid complex bodies (~3.3 yr) is 
close to a 7:2 mean motion resonance with Jupiter (3.39 yr), 
which has a shepherding effect, restricting the spread of 
dense swarms of meteoroids generated by disintegrations 
and so prolonging their lifetimes well beyond 20 yr (Asher 
& Clube 1993; Beech, Hargrove k. Brown 2004; Dubietis & 
Arlt 2007). From a study of historical meteor storms, and 
comparison with the cometary dust trails detected by IRAS, 
Kresak (1993) found that the survival times of the trails are 
typically ~60 yr. 

These considerations suggest that, in the course of dis- 
integration of the comet Encke progenitor over 20-30 kyr, 
one or more damaging encounters with the debris from a 
recent fragmentation is a reasonably probable event. 



Table 2. Orbital elements of Lake Tagish bolide (Brown et al. 
2000) compared to those of asteroid Cuno (Terentjeva & Bara- 
banov 2004) and Comet Encke. The elements are for epoch 2000. 



object 


a (AU) 


e 


i 


zo 


T{yr) 


Tagish Lake 


2.1 


0.57 


1.4 


160 


3.0 


(4183) Cuno 


2.0 


0.64 


6.8 


171 


2.8 


2P/Encke 


2.2 


0.85 


11.8 


161 


3.3 



4 SOME GEOPHYSICAL EXPECTATIONS 

We consider here some broadbrush expectations arising from 
the encounter of a swarm of cometary debris impinging on a 
hemisphere of the Earth. To match the observations at the 
Boundary, the swarm model must be able to reproduce the 
observations of extensive confiagration on at least a con- 
tinental scale, predict the high concentration of nanodia- 
monds, and account for the presence of apparently lunar 
material in the boundary layer. 

4.1 Thermal effects 

Wood subjected to a flux of radiant energy 
~4x lO^erg cm~^ s~^ will ignite within 1-20 seconds. 
While the whole-Earth encounter may take hours, the 
energy within any local horizon may be released overhead 
within a few seconds. Wildfires will thus be initiated over a 
land area say 5% of the area of the Earth (that of North 
America, for example) if the meteoroid swarm deposits ~1 
- 20x10^^ ergs during the encounter. At the extreme, we 
can assume an encounter speed d=30 kms~^ and 100% 
efficiency of conversion of kinetic to radiant energy in the 
upper atmosphere. Then the required mass of swarm to 
yield a continent- wide confiagration is ^2 - 40x10^^ gm. 
If the radiant energy is an order of magnitude down on 
this (say with the most of the energy going into bulk 
motion of air before disintegration, slowing the bolides to 
say 10 kms^^), then one requires the Earth to intercept 
a mass ~2 - 40x10^^ gm which is in the range of the 
masses discussed in Section 3. In reality, this energy will 
be distributed over a large number of discrete areas corre- 
sponding to Tunguska-like fireballs. There is therefore an 
expectation of conflagration over at least continental scales 
following such an encounter. A similar figure is reached 
if one distributes say N epicentres each of 2x10'' km^ 
of fire damage over 25x10^ km^ of continent: then one 
requires A'^ ~10* Tunguska-like impacts corresponding to 
an intercepted mass --^^SxlO^* gm. 

4.2 Nanodiamonds 

The Tagish Lake meteorite fell over Yukon Territory in 
British Columbia on 18th January 2000. The fragments are 
extremely fragile, low-density objects with a unique lithol- 
ogy, bearing some resemblance to the most primitive car- 
bonaceous chondrites types, the CI and CM chondrites; 
the meteorite is nevertheless quite distinct from either of 
them. It contains unprocessed pre-solar grains and the high- 
est abundance of nanodiamonds of any meteorite, ^^3650- 
4330 ppm (Grady, Verchovsky, Franchi et al. 2002). The 
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pre-entry orbital elements of the body are shown in Tabled 
along with those of asteroid Cuno and Comet Encke. The 
Tagish Lake fall has been identified with the /^-Orionid fire- 
ball stream which, in turn, has been associated with (4f83) 
Cuno (Terentjeva & Barabanov 2004), an asteroid belong- 
ing to the group of large NEOs which appear to be part 
of the Taurid Complex (Table 1). These generic relation- 
ships, along with the primitive nature of the Tagish Lake 
meteorite, suggest that it may be part of this Complex and 
that nanodiamonds are generally present in the progenitor 
comet. 

The concentration of nanodiamonds in the 12.9 ka sed- 
iments at multiple locations across North America ranges 
from ~0.01 to 3.7 ppm by weight (Kennett et al. 2009). Let 
the incoming cosmic material have fractional nanodiamond 
abundance /, falling on an area A of the Earth, and spread 
through a depth h of the Younger Dryas Boundary. Then to 
yield a concentration c throughout the boundary, the mass 
of swarm material intercepted by the Earth is 

M = Ahpc/f (1) 

where p is the bulk density of the material wherein the nan- 
odiamonds are found. If the bulk of the swarm material falls 
on say 5% of the Earth's surface {A = 2.5 x 10^^ cm^), then 
adopting / = 4 x 10"^ by weight from Grady et al. (2002), 
c = lO"'* from Kennett et al. (2009), and assuming the nan- 
odiamonds are concentrated in a layer of thickness h — 2 cm 
and bulk density p — 2.5 gmcm~^, one finds that intercept- 
ing a mass of 6 x 10^* gm can account for the observed 
nanodiamond concentration at 12.9 ka. This should proba- 
bly be increased by a factor of two or so because volatiles 
in the cometary material would disperse in the atmosphere, 
but it is still a modest requirement in the context of the 
breakup model. There is therefore no necessity for the nan- 
odiamonds to be formed in impacts, although of course this 
is by no means excluded. 

4.3 Lunar debris 

The Earth-Moon distance is ~60 Earth radii, comparable 
with the width of the meteoroid filaments being considered 
here (Fig. ^ ; thus any swarm which strikes the Earth may 
also strike the Moon, the bulk of this impacting mass be- 
ing in the few largest objects. The mass of lunar material 
thrown into space by impactors say > lO^'^ gm is uncertain 
but probably comparable to the mass of impacting material 
(Chyba, Owen & Ip 1994), and a proportion of this escaped 
material will strike the Earth. The dynamical evolution of 
lunar ejecta has been followed by Gladman, Burns, Duncan 
et al. (1995). The delivery efficiency p depends strongly on 
the location of impact, but characteristically p '^O. 05-0. 2 of 
the lunar material ejected will impact the Earth, the bulk of 
it within a few years of impact. Thus, if say the mass of the 
meteoroid swarm material striking the Earth is Me, then 
because of its smaller target area, the Moon will intercept 
0.05Af_B of swarm material, resulting in the ejection to space 
of a similar mass of lunar material and the fall of 0.05pM_B 
of Moon dust on to the Earth. With p = 0.1, one expects an 
admixture of lunar material at the wildfire horizon amount- 
ing to ~1% that of swarm material, to order of magnitude. 
This seems to be in the ballpark of the KREEP-like material 
discussed by Firestone (2009). 



5 SUMMARY AND CONCLUSIONS 

The proposition that an exceptionally large comet has been 
undergoing disintegration in the inner planetary system goes 
back over 40 years (Whipple 1967), and the evidence for 
the hypothesis has accumulated to the point where it now 
seems compelling. Radio and visual meteor data show that 
the zodiacal cloud is dominated by a broad stream of largely 
cometary material which incorporates an ancient, dispersed 
system of related meteor streams. Embedded within this sys- 
tem are significant numbers of large NEOs, including Comet 
Encke. Replenishment of the zodiacal cloud is sporadic, with 
the current cloud being substantially overmassive in relation 
to current sources. The system is most easily understood as 
due to the injection and continuing disintegration of a comet 
50-100 km in diameter. The fragmentation of comets is now 
recognized as a major route of their disintegration, and this 
is consistent with the numerous sub-streams and comoving 
asteroids observed within the Taurid Complex. The probable 
epoch of injection of this large comet, ~20-30 kyr ago, com- 
fortably straddles the 12.9 kyr date of the Younger Dryas 
Boundary. 

The hypothesis that terrestrial catastrophes may hap- 
pen on timescales ~0.1-1 Myr, due to the Earth running 
through swarms of debris from disintegrating large comets, 
is likewise not new (Clube & Napier 1984). However the 
accumulation of observations has allowed us to build an as- 
tronomical model, closely based on the contemporary en- 
vironment, which can plausibly yield the postulated YDB 
catastrophe. The interception of ~10^^ gm of material dur- 
ing the course of the disintegration is shown here to have 
been a reasonably probable event, capable of yielding de- 
struction on a continental scale. 

The object of this paper is not to claim that such an en- 
counter took place at 12,900 BP - that is a matter for Earth 
scientists - but to show that a convincing astronomical sce- 
nario can be constructed which seems to give a satisfac- 
tory match to the major geophysical features of the Younger 
Dryas Boundary data. 

If indeed the YDB event was an astronomical catastro- 
phe, its occurrence bears little relation to current impact 
hazard assessments derived from NEO surveys. 
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